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Preface 


The  purpose  of  this  thesis  was  to  determine  the 
feasibility  of  using  the  4x4  matrix  method  developed  by  Dwight 
Berreman  in  a  Spectroscopic  Ellipsometry  model  to  examine  the 
optical  properties  of  organic  and  polymer  thin  films.  When  it 
was  presented  as  a  possible  thesis  project,  it  sounded  too 
good  to  be  true.  Unfortunately,  that  turned  out  to  be  the 
case.  It  was  more  difficult  than  it  sounded,  much  like  the 
entire  AFIT  program.  My  getting  married  between  the  fourth 
and  fifth  quarter  didn't  make  things  any  easier  (at  school  or 
at  home) . 

There  was  more  coding  than  you  could  shake  a  stick  at  (in 
MATLAB,  no  less).  Considering  I  hadn't  written  any  code  since 
I  was  12  (and  that  was  in  BASIC),  it's  a  wonder  I  didn’t  go 
total  1 y  bald.  Fortunately,  my  wife  Diane  and  my  thesis 
advisor  Captain  James  D.  Targove  kept  me  from  going  insane. 
They  made  me  focus  on  one  problem  at  a  time.  And  for  that  I 
am  eternally  grateful. 

I  must  also  thank  the  Materials  Laboratory  of  the  Wright 
Research  and  Development  Center  for  providing  the  equipment 
and  samples  necessary  to  complete  this  thesis. 

Mark  L.  Sward,  Esq. 
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AFIT/GEP/ENP/ 90D-7 


Abstract 

This  thesis  demonstrated  the  feasibility  of  determining 
the  optical  properties  of  organic  and  polymer  thin  films 
through  the  use  of  Spectroscopic  Ellipsometry  (SE).  Tan  ¥  and 
cos  A  data  from  300-800  nanometers  (nm)  were  taken  with  a 
Rudolph  Research  s2000  spectroscopic  ellipsometer  on  four 
samples:  indium  tin  oxide  (ITO)  coated  glass;  five  layer 
poly-benzyl -L  glutamate  (PBLG)  organic  film  on  ITO  coated 
glass;  eight  layer  PBLG  film  on  ITO  coated  glass;  and  a 
thiophene  polymer  film  on  a  microscope  slide.  The  data  sets 
were  fit  to  a  choice  of  four  computer  models  based  on  a  paper 
written  by  Dwight  Berreman  in  1972.  The  four  models  were 
written  in  MATLAB  to  take  advantage  of  its  matrix  manipulative 
capabilities.  The  models  were:  a  single  layer  isotropic  film 
on  an  isotropic  substrate;  a  single  layer  anisotropic  film 
(with  random  orientation  of  optical  axes)  on  an  isotropic 
substrate;  two  isotropic  films  on  an  isotropic  substrate;  and 
two  anisotropic  films  (with  random  orientation  of  optical 
axes)  on  an  isotropic  substrate.  Using  only  tan  f  data  over 
a  restricted  wavelength  region,  all  four  data  sets  were  fit  to 
variances  of  0.01  or  less. 
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I .  Introduction 


Thin  films  (thicknesses  less  than  5000  A)  have  been 
recognized  and  used  for  well  over  a  century  now.  Whether 
naturally  occurring  or  man-made,  thin  films  affect  the 
electro-magnetic  properties  of  almost  all  materials  used  in 
today's  advanced  science.  Recently,  the  use  of  organic  thin 
films  (including  liquid  crystal  and  polymer  thin  films)  has 
been  considered  in  non-linear  signal  processing  as  switches, 
modulators,  and  wave-guides (1 ) .  In  order  to  predict  the 
properties  of  these  non-linear  devices,  one  must  first  have  a 
good  grasp  of  the  optical  properties  of  the  thin  films 
themselves . 

The  optical  properties  of  interest  are  usually  the 
complex  index  of  refraction  (N=n-ik)  or  the  complex 
dielectric  function  (e=ej  +  iej)  in  a  frequency  region  of 
interest.  An  important  way  of  obtaining  these  parameters  for 
thin  films  is  el  1 ipsometry .  Ellipsometry  is  a  non-destructive 
technique  that  measures  the  relative  change  in  polarization 
upon  reflection  from  a  sample.  The  reflectance  of  a  sample 
depends  on  many  parameters:  the  frequency  of  the  incident 
light;  the  angle  of  incidence;  the  number  of  films  on  the 
substrate;  the  thickness  of  the  films;  the  indices  of 
refraction  (or  dielectric  functions)  of  the  films  at  that 
frequency  and  that  angle  of  incidence;  whether  the  films  are 
anisotropic;  the  orientation  of  the  optical  axes  with  respect 
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to  the  film/substrate  surface;  the  optical  activity  of  the 
films;  surface  roughness;  whether  the  films  are  homogeneous; 
and  so  on. 

Once  the  reflectance  measurements  are  in  hand,  the 
optical  properties  and  thickness  of  any  thin  films  are 
determined  through  computer  modeling  of  the  sample.  The 
models  allow  for  a  certain  structure  of  the  film  and  include 
as  many  parameters  as  possible.  The  parameters  are  adjusted 
through  a  non-linear  least-squares  fitting  routine  until  the 
model  reflectance  matches  the  measured  reflectance  (within 
acceptable  tolerance). 

So,  in  order  to  understand  the  optical  properties  of 
organic  thin  films,  one  must  have  a  complete  description  of 
the  complex  indices  of  refraction  (or  dielectric  function) 
over  an  extended  frequency  region.  One  also  must  know  how  the 
films  are  oriented  with  respect  to  each  other,  or  the  sample 
surface,  in  order  to  obtain  the  desired  optical  response. 
Unfortunately,  the  ellipsometry  models  constructed  to  date  for 
organic  thin  films  have  not  taken  these  needs  into  account. 
Typically,  the  models  allow  for  only  one  frequency;  and 
therefore  only  one  value  of  the  index  of  refraction  (or 
dielectric  function)  is  calculated. 

This  thesis  was  undertaken  to  correct  some  of  these 
deficiencies.  Four  models  were  constructed  to  allow  for  the 
determination  of  the  complex  indices  of  refraction  across  an 
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optical  frequency  range.  These  models  were  successfully 
tested  against  published  data.  Finally,  the  models  were  used 
to  reduce  data  taken  from  multiple  layer  organic  films  (five 
and  eight  layers  of  poly  benzyl  L-glutamate  ( PBLG) ) .  These 
samples  were  mounted  on  an  indium  tin  oxide  coated  glass 
substrate.  The  models  were  also  used  on  a  film  of  thiophene 
mounted  on  a  microscope  slide. 

Although  multiple  wavelength  (frequency)  and  multiple 
angle  of  incidence  data  were  taken  and  successfully  reduced 
(as  mentioned  above),  the  models  used  were  hardly  complete. 
Three  of  the  model's  inherent  capabilities  (accounting  for 
optical  rotation,  paramagnetic  effects  and  random  orientation 
of  the  optical  axes)  were  not  taken  advantage  of.  Maybe  most 
importantly,  the  models'  accounted  for  neither  surface 
roughness  nor  inhomogeneous  indices  of  refraction. 
Nonetheless,  the  assumptions  made  in  modeling  the  organic  thin 
films  here  are  sound;  and  the  calculated  data  match  the 
measured  data  well. 
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I I .  Background 


El  1 ipsometry  measures  the  change  in  polarization  of  light 
upon  reflection  from  or  transmission  through  a  sample.  This 
thesis  dealt  exclusively  with  reflection  el lipsometry .  The 
change  in  polarization,  as  measured  by  an  el  1 ipsometer ,  is 
typically  characterized  by  the  two  parameters  tan  f  and  cos  A. 
Tan  ¥  represents  the  relative  amplitudes  of  the  perpendicular 
component  of  reflection  (Rs) ,  and  the  parallel  component  of 
reflection  (Rp).  Cos  A  is  the  difference  in  phase  between 
these  two  components.  The  defining  equation  is 

p  -  -  -j5e|  -  tan**14  (l) 


Note  that  the  reflection  coefficients  can  be  complex.  Good 
precision  is  expected,  since  tan  f  and  cos  A  are  parameters 
based  on  a  relative  measurement. 

Previous  el  1 ipsometric  studies  of  organic  thin  films  have 
concentrated  on  biological  samples  (2), (3).  Although 
interesting,  few  biological  samples  have  properties  analogous 
to  the  polymer  and  liquid  crystal  thin  films  needed  to  produce 
non-linear  optical  devices.  The  exceptions  are  biological 
samples  produced  as  Langmuir-Bl odgett  (LB)  films.  These  are 
films  in  which  the  individual  molecules  typically  attach 
themselves  to  the  substrate  at  one  end,  and  the  molecule 
itself  is  more  or  less  perpendicular  to  the  substrate  surface 
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(see  Figure  1 ) . 


Figure  1.  Typical  Molecule  Orientation  in  Langmuir-Blodgutt  Thin 
Films 


Biological  LB  films  do  resemble  the  polymer  LB  films 
necessary  for  non-linear  effects.  They  are  deposited  in  the 
same  fashion,  and  both  are  best  modeled  as  uniaxial  (two 
different  indices  of  refraction  in  different  directions)  thin 
films.  However,  previous  el lipsometric  studies  of  biological 
LB  films  (4)  restricted  themselves  to  null  el  1 ipsometry . 
Typical  null  el  1 ipsometry  has  a  polarizer,  compensator,  sample 
and  analyzer  (PCSA)  setup.  The  compensator  is  adjusted  such 
that  the  measured  reflectance  at  each  angle  of  incidence  and 
each  frequency  is  zero.  (For  a  complete  description  of  null 
el lipsometry ,  see  Azzam  and  Bashra  (5).)  Unfortunately, 
multiple  frequency  determinations  are  exceedingly  slow  in  null 
el lipsometry .  Therefore,  null  ellipsometry  lacks  any 
realistic  chance  of  yielding  the  necessarily  complete  optical 
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properties  of  polymer  or  liquid  crystal  thin  films.  A 
different  type  of  ellipsometry  is  required. 

A  Rotating  Analyzer/Polarizer  El  lipsometer  (RAE)  measures 
the  reflected  intensity  from  a  sample  and  Fourier  decomposes 
the  signal  into  a  calculated  tan  ¥  and  cos  A.  This  technique 
does  not  require  a  ''null"  determination.  By  taking  advantage 
of  the  computer  required  for  the  Fourier  decomposition,  one 
can  cover  a  large  frequency  range  in  a  short  period  of  time. 
When  using  RAE,  one  measurement  ii  obtained  at  each  frequency, 
with  two  resulting  data  points  (tan  ¥  and  cos  A).  RAE  is  the 
technique  required  to  completely  determine  optical  properties 
over  a  frequency  spectrum. 

Once  tan  ¥  and  cos  A  at  some  frequency  are  in  hand,  the 
optical  properties  and  thickness  of  any  thin  films  are 
determined  through  computer  modeling  of  the  sample.  The 
models  allow  for  a  certain  structure  of  the  film,  to  include 
as  many  parameters  as  possible.  In  RAE,  the  parameters  can 
include  indices  of  refraction,  orientation  of  optical  axes, 
optical  rotation,  paramagnetic  effects,  film  thicknesses  and 
surface  roughness. 

In  order  to  best  model  an  organic  thin  film,  this  thesis 
uses  the  matrix  technique  of  Berreman  (5).  This  models  the 
optical  properties  of  a  thin  film  as  a  6x6  matrix,  eventually 
yielding  the  4x4  reflection  matrix  for  the  entire  sample.  The 
model  assumes  the  geometry  of  Figure  2,  with  the  plane  of 
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incidence  the  x-z  plane. 


y  (out  _  x 

Film  (  s  ) 

Substrate  ^ 

Z 

Figure  2.  Assumed  Geometry  of  Sample 


The  derivation  of  the  optical  matrix  begins  by  casting 
Maxwell's  equations  in  a  6x6  matrix  (assuming  rationalized  MKS 
units  and  an  e1#t  time  dependence): 


0 


0 

0 


0 

_a 

dz 

_a 

dy 


0 

0 

0 

d_ 

dz 

0 

dx 


K 

Dx 

Dy 

-  io>  • 

Bz 

Bx 

By 

Bt 

(2) 


Substituting 
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0  0  0 

^  -  0  0  0 

0  0  0 


f  0  -±  ± 

dz  dy 

.  a  a 

“  &  0 

-  *  4  0 

dy  a x 


then  Eq.  (2)  can  be  rewritten  as 


cud Ez 

9  Hx 


or,  shortened,  as 


Q  •  G  -  to  •  C 


In  the  absence  of  nonlinear  optical  effects  and 
inhomogeneity ,  the  constitutive  relation  between  C  and  G  can 
generally  be  put  as: 
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£  -  MG 


(7) 


where 


(8) 


M  is  the  matrix  that  completely  characterizes  the  optical 
properties  of  the  thin  film.  u.  is  the  magnetic  permeability 

i 

tensor;  j>_  and  g_  are  optical  rotation  tensors;  and  e.  is  the 
dielectric  permittivity  tensor.  Allowing  for  an  anisotropic 
film  with  the  optical  axes  randomly  oriented  (after  rotation 
through  Euler  angles  a,  (3  and  y)  as  in  Figure  3,  e.  is  given 
by 

fi  -  -aL  A  (9) 

where  e3,  e2 ,  and  e3  are  the  principal  dielectric  constants 

ex  o  o' 

£  -  o  e,  o  (10) 

o  o  «3 
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and  the  coordinate  rotation  is  given  by  (6) 


A 


cosycoap-coaaainpainy  eoayain^+ooaacoaMiay  atnaainy 
-ainycoap-cosaainfcoay  -Binyalnft+coaacoaPcoBy  alnacoay 
ainaain0  -alnaooaP  eoaa 


(H) 


If  the  spatial  and  time  dependence  of  G  are  separated 

into 

a  -  e**  •  E  (12) 

then  Eq.  (6)  can  be  rewritten  as 

q  *  r  -  m  •  r  (i3) 


Look  again  at  Figure  2.  From  the  symmetry  of  the 
problem,  there  is  no  variation  in  the  y  direction  of  any  field 
component,  so 


a 


-  o 


(14) 


For  the  tangential  fields  to  match  across  the  boundary  at  z  =  0, 

-g  -  -l|»,  sin  e„  (15) 

where  <*  is  the  frequency  of  the  incident  light;  c  is  the  speed 
of  light  in  a  vacuum  (in  rationalized  MKS  units,  c  =  1);  Nq  is 
the  index  of  refraction  from  the  ambient  media;  0q  is  the 
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angle  of  incidence.  Let 


5  -  «  N0  sin60 


Then  Eq.  (4)  can  be  rewritten  as 


curl  - 


0  -i  0 

I  i  0  * 

.  0  -il  0. 


This  eliminates  Ex  and  Hx  from  Eq.  (13).  What  remains  is 


3  Hy 

JL  y  -  — i<i>  A  •  1 

dz  Ey  A  Ey 


X  -  -to A  •  X 


The  components  of  A  are  derived  from  M,  given  by 
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(20) 


An 

Aia 

An 

a14 

Aji 

A„ 

A„ 

Aj4 

Aji 

A,a 

Ajj 

AM 

A4i 

A«a 

A43 

A44 


and 


At  -  (M6lM36  -  M3XM„)/D 

A2  -  l(M„-n  )J#36  -M32M66]/D 

A3  -  (M€<M36  -  MUM")/D 

At  m  [M65M36  -  (Af35+t^)  Jfss] /D 

^5  -  (J#63M31  -  m33m61)/d 

At  -  [Mt3M32  -  (Mt2-i\)M33] /D  (21) 

^  -  (M<3j#34  -  m33m„)/d 

Ab  -  [(M35+ti)M„  -  M33M65]/D 

D  -  M33Um  -  M3tM63 

n  -  —  sin0o 

(4  C  0 


With  A  known  from  M,  »,  Ng,  and  6g,  one  can  predict  how  the 
electric  and  magnetic  waves  will  travel  in  the  film. 

In  the  general  case,  M  is  an  arbitrary  function  of  z  and 
Eq.  (19)  cannot  be  solved  analytically.  However,  assuming  the 
film  is  homogeneous,  M  is  independent  of  z.  Then  Eq.  (19)  can 
be  integrated  directly  to  give 
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I(z  +  h)  -  Uh)  -I(z) 


(22) 


where 

L(h)  -  (23) 

is  the  4x4  "layer"  matrix  for  a  thin  film  of  thickness  h.  Eq. 
(22)  gives  the  ability  to  model  the  reflection  from  a  sample 
as  a  matrix  and  include  it  in  el lipsometric  calculations  (as 
will  be  shown  later). 
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If  more  film  layers  are  present:  add  more  "layer" 
matrices;  multiply  the  "layer"  matrices  together;  and  use  the 
resulting  matrix  as  described  below.  Specifically,  if  the 
sample  is  to  be  modeled  as  in  Figure  4,  the  overall  layer 
matrix  would  be  given  by  L  =  LXhj^LXhj ) . 


Film  1  ( hi ) 

Film  2  ( h 2 ) 


Substrate 


Figure  4.  Modeling  with  Two  Thin  Films  on  a  Substrate 

As  you  can  see,  Berreman's  technique  allows  for  films  to 
be  modeled  as  anisotropic,  optically  rotating,  paramagnetic, 
and  homogeneous.  All  of  the  models  used  in  this  thesis 
assumed  the  following: 

1)  Films  were  non-magnet ic  (y=l);  smooth;  homogeneous; 

and  have  no  optical  rotation. 

2)  The  ambient  refractive  index  (Nfl)  was  1.000. 

3)  All  films  were  on  isotropic  substrates. 
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Ill .  Computational  Methods 

It  is  still  unclear  how  one  gets  from  the  layer  matrix 
calculation  (Eq.  (23))  to  the  el  1 ipsometric  parameters  tan  f 
and  cos  A.  A  diagram  of  the  el lipsometric  configuration  used 
in  this  experiment  is  given  in  Figure  5: 


First,  model  the  light  through  each  optical  device  as  a 
4x4  matrix.  Assume  the  light  striking  the  fixed  polarizer  xs 
entirely  unpolarized  (untrue;  but  it  will  be  shown  the  fixed 
polarizer  completely  polarizes  the  light).  The  Stokes  vector 
describing  the  light  is  given  by 


16 


(24) 


s.  - 


T 

o 

0 

p. 


After  passing  through  the  fixed  polarizer,  the  Stokes 
vector  is  now 


5-4 


10  10 

1 

1 

0  0  0  0 

0 

1 

0 

10  10 

0 

2 

1 

0  0  0  0. 

0. 

0 

(25) 


After  passing  through  the  second  polarizer  (rotating  at 
an  angular  velocity  of  *),  the  resulting  Stokes  vector  becomes 
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1 

1 

cos20 

sin2© 

0 

0 

.  1 

COS2GJ 

COSz2C3 

sin2cscos2CJ 

0 

1 

2 

sin2Q 

sin2C9cos2C3 

sina2t0O 

0 

p, 

0 

0 

0 

0. 

(26) 
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1  +  sin2c) 

2  8in2C9  COS2C3  +  COS2C3 
4  sina2©  +  sin20 
0 


(27) 


The  next  component  in  the  optical  path  is  the  sample.  If 
the  sample  surface  is  assumed  to  be  homogeneous,  it  is 
characterized  by  el  1  ipsometric  parameters  tan  ¥  and  cos  A,  its 
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Mueller  matrix  is  given  by  (5): 


£  - 


(tana$+l) 

(tana^-l) 

0 

0 


(tan2f-l)  0  0 

(tan2t+l)  0  0 

0  2cosAtanf  2sinAtanf 
0  -2sinAtanf  2cosAtanf. 


(28) 


So,  the  Stokes  vector  after  reflection  from  the  sample 
surface  is  given  by 


& 


(l+sin2ttf)  (tan2$+l)  +  (cos2eJ+sin2t8cos2«J)  (tana$-l) 
l  (l+sin2ta)  (tan2<r-l)  +  (cos2CJ+sin2taco82ta)  (tan2%+l) 

4  (sin22ta+sin2a)  (2cosAtan$) 

(sina2GJ+sin2c3)  (-2sinAtanf) 

(29) 


The  light  passes  through  the  fixed  analyzer,  producing 
the  Stokes  vector  that  reaches  the  photomultiplier  tube  ( PMT ) . 
The  PMT  measures  the  intensity,  given  by  S(l): 


—  [(l+sin2W)  (tan2^+l) 

8 

+  (cos20 +sin20cos2f9)  (tana^-l) 
+  (2cosAtant)  (sina2t3+sin2©)  ] 


(30) 


Substituting 


sin2Gcos2C3 

sina2Q 


~sin4C3 

2 

~  -  -i-cos4Q 
2  2 


(31) 
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What  is  needed  is  the  4x4  Mueller  matrix  for  the  sample,  M. 
It  is  not  the  same  as  the  layer  matrix,  L.  The  technique  is 
to  go  from  the  4x4  layer  matrix  L,  to  the  2x2  sample 
reflection  matrix  R,  where 


19 


R  - 


K  *« 


to  the  4x4  sample  Mueller  matrix  M.  For  any  layered  structure 
on  an  isotropic  substrate  (4),  R  is  given  by 


R  - 


where 


N0  (^jA^-J^jCosej) 
N0  (L^-^cosdj) 
(1^3^-12300863)  + 
(133^,-13300863)  - 
(L~32Neco8B2-L42) 
N0  (L33NMcoa&2  -L42 ) 
<133^00863-143)  + 
<133^^00863—143)  — 


+  cos60 

-  ro860 
tf0cos60 
N0cob60 
+  cos60 

-  cos60 
AT0cos80 
N0cobBq 


(L^-^cobB^ 
<li  i^.-l2ic°8®2) 
<1241^-13400862) 
<  863  ) 

(IjjNjCOSea-^) 
<132^00863-141) 

(L34,Njco8B2-L44) 
<134^00863-  Au> 


Now,  the  4x4  sample  Mueller  matrix  is  given  by  (4): 


£  -  (*®&1)  *4^; 


where 


1  - 


10  0  1' 
10  0-1 
0  110 
0H0 
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R  the  2x2  sample  reflection  matrix;  R*.  is  the  complex 
conjugate  of  R;  and  R  •  r1  is  a  direct  product,  given  by 
(8): 


&1 


^12*^ 

Rix'Rl  R22mBl 


(41) 


So,  the  Stokes  vector  after  reflection  from  the  sample 
surface  is  given  by 

+fifia  (8in20co82W+cos20)  (sina2048in2® ) 

1  St2l  (l>sin20)  +&22  ( Bin20cos20+cos20 )  *tt33  (sina2tD+Bin2©) 

*  «33  (1+Bin20)  ^32  (8in20cos20-f cos20)  +H,j  (sina20+ sin20)  I 
H41  (l+sin2«)  +H4a  (8in20cos204cos20)  +Uf43  (8ina2«+sin2«) 

The  light  passes  through  the  fixed  analyzer,  producing 
the  Stokes  vector  that  reaches  the  photomultiplier  tube  (PMT) : 

lMxl+Mn)  (i+sin2w)  *  (sln2«»coa2v«co«2e)  +  (tu«Sul  (»in*2»M«in2«) 

0 

(8ix*tt,x)  <l+*in2»)  (*1ji2»co»2»+co«2«)  ♦  (Mx,*8„)  (■in328+»in2«) 

0 


Using  Eq.(31),  and  recalling  that  S(l)  is  rewritten  as  a 
truncated  Fourier  series. 
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+  ^  COB2a 

O 

5(1)  -  -  ~  (H^+lljj)  C0S4O  (44) 

♦  -|  ain2« 

+  (H12+#3a)  sin4ta 

Comparison  of  Eq.  (44)  with  Eq.  (33),  we  see  that 

“  "g  ^ll+^31+'g  (^j+AIjj)  ] 

"  *g  (^13+^33) 
fl4  “  ^13+^33^ 

^2  "  "g  (^ll+^31+^13+^33 ) 

^4  "  (&i2+&n) 


Substituting  the  Fourier  coefficients  of  Eq.  (45)  into  Eqs. 
(34)  and  (35),  the  calculated  tan  f  and  cos  A  are  found. 

Let  the  sample  be  measured  over  some  wavelength  region, 
giving  a  set  of  tan  ¥  and  cos  A  for  each  wavelength.  Assume 
that  a  model  with  m  parameters  (indices  of  refraction, 
thickness,  etc)  was  chosen.  How  is  it  known  when  the  best 
estimate  of  the  model  parameters  is  obtained?  More 
importantly,  how  is  an  initial  guess  at  each  of  the  m  model 
parameters  improved  to  get  to  the  best  estimate? 

The  usual  method  of  minimizing  the  weighted  sum  of  the 
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squared  deviations  (least-squares  fit)  is  followed.  Assume 
all  measurement  errors  follow  a  normal  (Gaussian) 
distribution.  Then  define  a  mean  squared  error,  %K  as: 


Xa 


1 

N- 1 


(46) 


where  N  is  the  number  of  data  points;  y^ex  is  the 
experimentally  determined  value  of  either  tan  ¥  or  cos  A;  y^c*^c 
is  the  calculated  value  (from  the  model)  of  tan  ¥  or  cos  A; 
and  is  the  standard  deviation  of  the  measured  value. 

According  to  the  method  of  least-squares,  the  optimum 
values  of  the  m  model  parameters  3  =  a^,  a2,  a3,  ...  ,  a^,  ... 
at  are  obtained  by  minimizing  x2  with  respect  to  each  of  the 
parameters  simultaneously 


X2  must  be  considered  a  continuous  function  of  the  m 
model  parameters  describing  a  hypersurface  in  m-dimensional 
space.  This  space  must  be  searched  for  the  appropriate 
minimum  value  of  y} .  One  the  difficulties  of  such  a  search 
is  there  may  be  more  than  one  local  minimum  for  y}  within  a 
reasonable  range  of  values  for  the  3. 

There  are  many  methods  of  searching  the  parameter  space. 
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Marquardt  derived  a  method  based  on  a  suggestion  by  Levenberg 
(9).  The  Marquardt-Levenberg  method  has  become  the  standard 
non-linear,  least  squares  fitting  routine.  It  varies  between 
two  different  methods:  "steepest  descent"  and  "inverse 
Hessian” . 

The  inverse  Hessian  method  expands  the  "fitting  function" 
y  as  a  function  of  the  m  model  parameters  3.  Suppose  y  is 
expanded  to  second  order  in  a  Taylor's  series  as  a  function  of 
a : 


(48) 


If  the  derivatives  of  x2  with  respect  to  the  model  parameter 
increments  63  are  set  equal  to  0,  and  neglecting  all  terms  of 
order  one  or  higher  in  6aj,  the  resulting  equations  are 
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where  (3  is  called  the  gradient  vector  and  a  is  called  the 
curvature  matrix.  Eq.  (48)  is  equivalent  to  approximating  the 
X2  hypersurface  with  a  parabolic  surface. 

The  problem  with  the  inverse  Hessian  method  (named 
because  the  matrix  a  is  a  Hessian  matrix)  is  calculating  the 
step  sizes  63.  The  step  sizes  must  be  large  enough  to  prevent 
round-off  error  in  the  computation,  and  small  enough  to 
furnish  reasonable  results  near  the  minimum  where  the 
derivatives  may  be  changing  rapidly  with  the  parameters. 

If  the  model  parameters  are  close  enough  to  the  minimum 
to  invoke  parabolic  approximation,  then  Eq.  (50)  can  be  solved 
directly  to  yield  the  parameter  step  sizes  63  such  that  x 
should  be  minimized  by  3+63.  If  the  starting  point  is  close 
enough  to  the  minimum  that  higher-order  terms  in  the  expansion 
can  be  neglected,  inverse  Hessian  is  an  accurate  and  precise 
method.  But  if  the  starting  point  is  not  close  enough  to  the 
minimum,  the  parabolic  approximation  is  invalid  and  the 
results  are  in  error. 

In  contrast,  the  steepest  descent  method  is  ideally 
suited  for  approaching  the  minimum  in  y}  from  far  away,  but 
does  not  converge  rapidly  when  in  the  immediate  neighborhood. 
In  this  method,  all  of  the  parameters  3  are  incremented 
simultaneously,  with  the  relative  magnitudes  adjusted  so  that 
the  resultant  travel  in  parameter  space  is  along  the  gradient 
(the  direction  of  maximum  variation)  of  x  • 
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The  gradient  of  X2,  Vx2,  is  given  by 


VX2 


(51) 


In  order  to  determine  the  gradient,  the  variation  of  x2  in  the 
neighborhood  of  the  starting  point  is  sampled  independently 
for  each  parameter  to  yield  an  approximate  value  for  the  first 
derivative 

W),  -  %  -  lis>  (52) 


Note  that  the  dimensions  (and  units)  of  all  of  the  model 
parameters  are  not  the  same.  This  is  solved  by  letting 


(53) 


Define  a  dimensionless  gradient  y 


(54) 


The  method  of  steepest  descent  follows  the  opposite  direction 
of  the  dimensionless  gradient  y 


baj  -  -y j  A aj  (55) 

Note  the  similarities  between  Eq.  (49)  and  Eq.  (54). 
Here  is  where  the  two  methods  are  combined  to  get  the  best 
features  of  both.  If  the  diagonal  of  the  curvature  matrix  a 
is  increased  by  a  factor  of  A,  Eq.  (49)  becomes 

fj  -  £  h  •*  (56) 


where 


Oj(l+X) 

** 


(j*k) 


(57) 


If  X  is  small,  Eq.  (56)  acts  like  an  inverse  Hessian.  If 
X  is  large,  the  diagonal  terms  of  the  curvature  matrix 
dominate  and  the  matrix  equations  become  gradient  equations 
(i.e.,  steepest  descent).  The  solution  for  the  step  sizes  in 
either  case  is  given  by 


where  a'-1  is  the  inverse  of  the  augmented  curvature  matrix. 
The  inverse  of  the  a  matrix  is  called  the  error  matrix.  It  is 
used  to  calculate  the  90%  confidence  limits  of  the  model 
parameters . 

Thus,  given  an  initial  set  of  m  parameters,  3, 

1)  Compute  X2(2); 

2)  Start  with  X  =  0.001; 

3)  Compute  x2(^+63); 

4)  If  x*(Z+6a)  >  x2(3),  increase  X  by  a  factor  of  10  and 

redo  3); 

5)  If  x*(^+63)  <  X^(3),  decrease  X  by  a  factor  of  10; 

update  the  values  of  3  =  3+63;  and  redo  3). 

Once  a  minimum  in  x2  is  reached,  the  set  3  contains  the 
parameters  that  best  fit  the  model . 

How  is  one  certain  the  set  of  3  parameters  are  truly  the 
best  fit?  The  error  matrix  allows  one  to  calculate  the 
standard  error  and  90%  confidence  intervals  for  each 
parameter.  The  standard  error  of  each  model  parameter  is 
given  by: 


and  the  90%  confidence  interval  for  each  model  parameter 
varied  one  at  a  time  is  given  by  (10) 
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ft,  -  o,  JT7JT  -  1.65  o, 


(60) 


Thus,  the  signs  of  a  good  fit  are: 

1)  small  x*  («1); 

2)  reasonable  90%  confidence  intervals  (i.e.,  not  tf  =  100  ± 
75  A,  or  nj  =  1.51  ±  1.25). 

However,  there  is  an  inherent  problem  with  SE  data  taken 
at  only  one  angle  of  incidence.  Assume  a  model  has  five 
parameters.  At  any  one  wavelength,  only  two  data  points  (tan 
V  and  cos  A)  are  obtained  for  each  measurement. 
Unfortunately,  the  laws  of  statistics  dictate  that  for  a  good 
fit  to  m  model  parameters,  at  least  (m+2)  data  points  must 
be  used.  Therefore,  multiple  angle  of  incidence  (MAI)  data  is 
needed  at  each  wavelength  in  order  to  obtain  a  good  fit. 

■ 

The  models  used  in  this  thesis  were  written  in  MATLAB  to 
take  advantage  of  its  powerful  and  straight-forward  matrix 
capabilities.  A  simplified  flow  chart  of  the  models  is  given 
in  Figure  6. 

Four  models  were  constructed,  with  the  following 
parameters : 

1)  One  isotropic,  absorbing  film  (n^(A),  kf(A),  t{ )  on  an 
isotropic,  absorbing  substrate  (n,(A),  ks(A)); 

2)  One  anisotropic,  absorbing  film  (n^A),  kj(A),  n^A),  ky(A), 

n,(A),  kt(A),  Of,  (3{ ,  tf)  on  an  isotropic,  absorbing 

substrate  (ns(A),  k,(A)); 
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input  data  file  containing 
6  ,  aaaatzed  tan  f ,  cea  A 


(rixtt  Ctii) 
fiat  X-  #.m 

Sat  jatiiti  foi  aodal  paiaaetais,  a  | 

Caaputa  Xta).  ,  ^(a) 

fiat  which  pazaaataia  aza  to  ba  railed  j 

JCoapute  a  +  <fa 

I 

I  Ceipita  <Ta) ,  Jjl;i»<fa),  .£>U*cfa) 

I 


t 


I 


I 


la  VltlAICl  >  • .  •  1  at  X<  lxie  ? 


Taa  ^ 


la 


*  -*■ 


"X  •  i o •  "X  .» 

*StMj 


\e*  *  1/10  ' 

Stoza  (a*  <Ta ) ,  ^  it*  (f*) ,  £*) ,  (a  ♦  cfa) 

Calculate  confidence  intertill 


Stoza  final  paranetezi,  tot  coafldanca  iatarrala 
Figure  6.  MATLAB  SB  Model  Flow  Chart 
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3)  Two  isotropic,  absorbing  films  (nfj(A),  k^(A),  tfl ,  nf2(A), 
k{j(X),  tjj)  on  an  isotropic,  absorbing  substrate  (ns(A), 
ks(A)); 

4)  Two  anisotropic,  absorbing  films  (n^A),  kxj(A),  n^(A), 
kjj  ( A ) ,  R||  ( A ) ,  kjj  ( A ) ,  oij ,  (3^  t  Yj  <  tj|  , n^  ( A ) ,  k^2  ( A ) ,  n^ 1 
kj2 ( A ) ,  nl2 ( A ) ,  kl2(A),  a2,  '  *2'  fcf2>  on  an  isotropic, 
absorbing  substrate  (nj(A),  ks(A)). 

The  basic  models  of  the  f ilm(s)/substrate  system  were 
checked  against  a  variety  of  published  data 
(11 ) , ( 12) , (13) , (14 ) , ( 15) , (16) .  The  results 

matched  very  well  and  validated  the  basic  code. 
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IV.  Experimental  Information 


All  of  the  data  used  in  this  thesis  was  taken  with  a 
Rudolph  Research  Corporation  s2000  Spectroscopic  Ellipsometer 
(7).  A  diagram  of  the  machine  layout  was  shown  in  Figure  5. 
The  light  source  used  in  the  s2000  was  a  75  watt,  high- 
pressure  Xenon  lamp.  The  fixed  polarizer  was  an  UV  grade 
Glan-Thompson  prism,  set  (for  all  thesis  data  taken)  at  45.00 
±  0.01*.  This  polarizer  converted  the  collimated  light  into 
linearly  polarized  light,  with  equal  parallel  (p)  and 
perpendicular  (s)  components. 

The  rotating  polarizer  was  also  an  UV  grade  Glan-Thompson 
prism.  It  converted  the  linearly  polarized  light  into  a  beam 
whose  intensity  and  polarization  direction  were  modulated  by 
the  rotation  of  the  prism.  The  fixed  analyzer  was  also  an  UV 
grade  Glan-Thompson  prism  set  at  45.00  ±  0.01*.  This 
converted  the  reflected  beam  into  an  output  beam  modulated  in 
intensity  only. 

The  monochromator  used  in  the  s2000  was  a  dual 
monochromator.  The  PMT  was  a  Hammamatsu  R928.  A 
preamp/digitizer  converted  the  current  output  from  the  PMT  to 
a  voltage  signal  and  stored  the  intensity  as  a  12  bit  word. 

The  s2000  measured  the  el lipsometric  parameters  tan  ¥  and 
cos  A  as  follows: 

1)  Set  the  monochromator  to  0; 

2)  Drove  the  monochromator  to  the  starting  wavelength; 
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3)  Adjusted  the  PMT  voltage  to  get  a  signal  between  5 

and  7.5  volts; 

4)  Averaged  measurements  of  tan  V  and  cos  A  until 

-  standard  deviation  less  than  0.001,  or 

-  100  measurements  were  made; 

5)  Stored  the  wavelength,  tan  V  and  cos  A  data;  and 

continued  to  the  next  wavelength. 

Four  samples  were  analyzed  for  this  thesis.  Sample  1  was 
a  indium  tin  oxide  (ITO)  coated  glass  substrate.  Sample  2  was 
five  layers  of  PBLG  on  a  substrate  of  ITO  coated  glass. 
Sample  3  was  eight  layers  of  PBLG  on  a  substrate  cf  ITO  coated 
glass.  All  of  the  substrates  for  these  three  samples  had  the 
same  dimensions:  10  millimeters  (mm)  in  length,  6  mm  in  width 
and  2  mm  in  thickness. 

This  substrate  thickness  was  too  small  for  the  s2000  to 
resolve  the  differences  between  the  front  and  back  surface 
reflections.  Since  directly  mounting  the  samples  on  the  s2000 
vacuum  chuck  didn't  solve  the  problem,  a  different  approach 
was  tried.  The  substrates  were  affixed  to  a  fused  silica 
(SiOj)  block  with  index  matched  NYOGEL*.  This  allowed  the 
front  and  back  surface  reflections  to  be  resolved.  A  test 
with  a  similar  thickness  of  BK7  glass  mounted  with  NYOGEL  to 
the  fused  silica  block  validated  the  technique. 

Sample  4  was  a  film  of  thiophene  mounted  on  a  microscope 
slide.  The  back  side  of  the  slide  was  spray  painted  black  in 
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order  to  eliminate  the  back  surface  reflections.  The  slide 
was  mounted  directly  onto  the  vacuum  chuck  and  the  film 
analyzed. 

Each  of  the  four  samples  were  analyzed  from  300  to  800 
nanometers  (nm)  at  every  5  nm.  Measurements  were  taken  across 
the  spectrum  at  60.00*,  65.00*,  70.00*,  75.00*  and  80.00*  for 
all  four  samples.  Plots  of  tan  ¥  and  cos  A  for  all  four 
samples  at  all  five  angles  are  given  in  Appendices  A-D. 


34 


V.  Results  and  Discussion 

Additional  data  was  taken  to  determine  the  repeatability 
of  the  measurements.  This  would  also  gave  a  range  to  which 
the  data  could  be  fit.  There  is  no  point  to  fit  the  data  to 
a  variance  of  0.0001  when  the  data  is  only  repeatable  to  0.01. 

Tw>.  sets  of  data  were  taken  on  the  5  layer  PBLG  film,  one 
immeuiately  after  the  other.  The  differences  between  the  two 
tan  f  and  cos  A,  every  5  nm  from  300-400  nm,  are  tabulated  in 
Table  1.  Note  that  the  average  change  in  tan  f  measurements 
was  0.0021. 

An  example  of  two  measurements  made  some  time  apart  on 
different  spots  is  given  in  Table  2.  Again,  two  data  sets 
were  taken  of  the  5  layer  PBLG  film  system.  They  were  taken 
apprc  :imately  30  minutes  apart,  with  the  substrate  in  the  same 
basic  orientation.  Again,  the  tabulated  data  are  the 
differences  between  the  two  tan  f  and  cos  A.  The  average 
change  in  tan  Y  measurements  was  about  0.0051. 

This  would  seem  to  set  a  convergence  criteria  of  about 
°.01.  In  other  words,  when  the  variance  is  less  than  0.01, 
then  we  can  assume  the  fit  is  as  good  as  another  data  set 
taken  on  the  same  sample. 
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TABLE  1 


MEASUREMENT  TO  MEASUREMENT  CHANGES  (5  LAYER  PBLG) 


cos  A) 

300 

0.0054 

0.0050 

305 

0.0039 

0.0021 

310 

0.0001 

0.0019 

315 

0.0001 

0.0023 

320 

0.0033 

0.0012 

325 

0.0051 

0.0010 

330 

0.0019 

0.0016 

335 

0.0056 

0.0018 

340 

0.0130 

0.0055 

345 

0.0651 

0.0021 

350 

0.0468 

0.0043 

355 

0.0951 

0.0015 

360 

0.0488 

0.0023 

365 

0.0257 

0.0007 

370 

0.0240 

0.0029 

375 

0.0158 

0.0015 

380 

0.0286 

0.0013 

385 

0.0127 

0.0019 

390 

0.0029 

0.0006 

395 

0.0140 

0.0014 

400 

0.0050 

0.0008 

TABLE  2 


TIME  TO  TIME  (OR  SPOT  TO  SPOT)  CHANGES  (5  LAYER  PBLG) 


300 

0.0134 

0.0046 

305 

0.0082 

0.0045 

310 

0.0196 

0.0078 

315 

0.0354 

0.0063 

320 

0.0390 

0.0055 

325 

0.0442 

0.0027 

330 

0.0388 

0.0056 

335 

0.0613 

0.0047 

340 

0.0590 

0.0034 

345 

0.0735 

0.0044 

350 

0.0684 

0.0020 

355 

0.0501 

0.0038 

360 

0.0707 

0.0046 

365 

0.0440 

0.0050 

370 

0.1176 

0.0044 

375 

0.0394 

0.0060 

380 

0.0415 

0.0031 

385 

0.0504 

0.0042 

390 

0.0174 

0.0083 

395 

0.0198 

0.0067 

400 

0.0191 

0.0095 

Unfortunately,  this  also  puts  a  restriction  on  the 
quality  of  the  fit  [see  Eq.  (59)].  For  example,  say  a  film  had 
a  small  degree  of  anisotropy.  Then  let  the  difference  in 
variance  between  data  fitted  with  an  isotropic  model  and  data 
fitted  with  the  correct  degree  of  anisotropy  be  smaller  than 
0.01.  The  "true  fit"  would  be  impossible  to  determine 
accurately  given  the  quality  of  data  described  above. 

However,  it  appears  from  the  data  in  this  thesis  that 
this  might  not  be  much  of  a  problem.  Note  that  the  PBLG  film 
layers  had  small  degrees  of  anisotropy  (0.03).  The  variances 
in  the  fits  to  these  data  sets  were  acceptable  (less  than 
0.01)  only  when  this  anisotropy  was  accounted  for. 

Another  restriction  placed  on  fitting  has  to  do  with 
dispersion.  Normally  in  the  visible  region  of  the  electro¬ 
magnetic  spectrum,  both  the  real  and  imaginary  parts  of  the 
index  of  refraction  for  organic  films  decrease  as  X  increases. 
However,  the  models  were  not  able  to  pull  a  dispersion  profile 
out  of  a  flat  guess  for  samples  that  were  known  to  be 
dispersive.  The  models,  did  however,  improve  a  dispersive 
guess . 

The  sample  measured  in  this  thesis  known  to  be  dispersive 
was  the  ITO  coated  glass.  Based  on  a  priori  knowledge  of  the 
ITO  sample,  it  was  determined  that  the  best  fit  would  come 
from  the  single  isotropic  film  model.  Examination  of  the  ITO 
data  showed  at  wavelengths  greater  than  370  nm,  cos  A  and  tan 
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f  results  were  not  strongly  influenced  by  the  film.  Use  of 
data  above  this  wavelength  would  be  worthless.  Hence  all  of 
the  data  used  for  the  determination  of  the  optical  properties 
of  the  ITO  films  (Samples  1-3)  was  from  305  to  370  nm. 

Additionally,  since  large  changes  in  phase  produce  small 
changes  in  cos  A  when  cos  A  is  greater  than  I  0.7071  ,  this  data 
is  not  normally  used  for  accurate  fitting.  For  some  angles  of 
incidence,  all  of  the  cos  A  data  was  greater  than  I  0.7071  . 
Therefore,  it  was  decided  to  use  only  the  tan  f  data  in 
reducing  the  data  for  all  four  samples. 

In  fine  tuning  the  guesses  made  to  fit  data,  some 
observations  were  made.  If  it  is  assumed  all  other  parameters 
are  fixed: 

1)  as  the  real  part  of  the  index  of  refraction  (n^)  went 
up,  tan  T  went  up  and  cos  A  went  down,  and  vice  versa; 

2)  as  the  complex  part  of  the  index  of  refraction  (kj) 
went  up,  tan  T  went  down  and  cos  A  went  up,  and  vice 
versa; 

3)  as  the  film  thickness  (tf)  went  up,  tan  T  went  down 
and  cos  A  went  up,  and  vice  versa; 

4)  for  all  of  the  samples  used  in  this  thesis,  the 
fitted  data  was  most  sensitive  to  film  thickness. 

Since  it  was  discovered  that  the  models  would  not  move 
automatically  to  the  "true"  minimum  no  matter  what  the  initial 
guess,  a  search  strategy  was  devised  taking  items  1-4  into 
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account . 


First,  since  the  models  were  most  sensitive  to  film 
thickness,  the  complex  part  of  the  index  of  refraction  would 
initially  be  set  to  zero  and  the  real  part  of  the  index  of 
refraction  fixed  at  the  same  "arbitrary"  value  for  all 
wavelengths.  Arbitrary  meant  within  10%  of  the  expected  value 
of  the  real  part  of  the  index  of  refraction.  Then  the  fit 
would  be  checked  at  different  film  thickness  near  the  expected 
value. 

Once  the  best  fit  with  a  given  film  thickness  was 
determined,  the  thickness  was  fixed  and  the  real  part  of  the 
index  of  refraction  was  varied.  The  fits  usually  changed  very 
little  with  a  "flat"  guess  at  index  of  refraction  and  fixed 
thickness.  Once  the  best  fit  with  a  flat  index  profile  and 
film  thickness  was  set,  the  experimental  data  was  examined  to 
see  if  absorption  appeared  to  be  present.  If  so,  then  the 
complex  part  of  the  index  of  refraction  was  set  at  some  small 
non-zero  value  and  allowed  to  vary. 

If  the  variance  was  still  not  less  than  0.01  after  all  of 
this,  a  choice  was  made.  In  the  case  of  films  modeled  as 
isotropic  (i.e.,  ITO) ,  some  dispersion  was  added  to  the 
initial  guess  at  the  indices  of  refraction  (both  real  and 
imaginary  parts  if  necessary).  In  the  case  the  films  were  to 
be  modeled  as  anisotropic  (i.e.,  PBLG  5  and  8  layer  films  on 
ITO),  then  the  values  from  the  isotropic  models  would  be  put 


40 


into  the  anisotropic  model . 

Dispersive  guesses  were  made  for  both  the  real  and 
imaginary  parts  of  the  index  of  refraction  for  the  ITO  coated 
glass.  The  final  fitted  parameters  produced  a  variance  of 
0.0093.  The  film  thickness  was  calculated  to  be  279  ±  3  A. 
The  measured  and  calculated  data  from  305-370  nm  at  all  angles 
for  the  ITO  coated  glass  are  plotted  in  Appendix  E. 
Additionally,  the  real  and  imaginary  parts  of  the  index  of 
refraction,  along  with  the  corresponding  90%  confidence 
intervals,  are  plotted  versus  wavelength  for  the  ITO  film  in 
Appendix  F. 

The  five  and  eight  layer  PBLG  films  were  fit  with  the  two 
layer  anisotropic  film  model.  They  were  both  tried  with  the 
two  layer  isotropic  film  model,  but  the  variances  in  the  fits 
with  the  tan  T  data  was  never  less  than  0.01.  The  smallest 
variance  with  the  two  layer  isotropic  film  model  was  0.014. 
The  final  variance  in  the  fit  for  the  five  layer  PBLG  film  was 
0.0075,  with  flat  indices  of  Nx  =  1.65  -  0.20*i,  Nj  ~  Nx  =  1.62 
-  0.10*i,  and  film  thickness  of  Tj  =  164  ±  2  A.  For  the  eight 
layer  PBLG  film,  the  final  variance  was  0.0056,  with  flat 
indices  of  Nx  =  1.65  -  0.18*i,  Nj  =  Nx  =  1.62  -  0.09*i,  and 
film  thickness  of  Tf  =  216  ±  2  A.  The  measured  and  calculated 
data  from  305-370  nm  at  all  angles  of  incidence  for  the  five 
layer  PBLG  film  are  plotted  in  Appendix  G.  The  measured  and 
calculated  data  from  305-370  nm  at  all  angles  of  incidence  for 
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the  eight  layer  PBLG  film  are  plotted  in  Appendix  H. 

The  thiophene  film  was  fit  with  the  single  isotropic  film 
model.  In  this  case,  neither  dispersion  nor  absorption  was 
needed  to  fit  the  data  to  within  0.0075.  Here,  the  tan  f  data 
from  500-800  nm  was  used.  The  data  taken  at  60*  angle  of 
incidence  was  not  included  in  the  fit  because  it  appeared  the 
back  surface  reflection  was  not  removed  by  the  black  paint  at 
that  angle.  The  film  thickness  was  calculated  to  be  3060  ±  7 
A,  with  an  uniform  index  of  refraction  from  500-800  nm  of  Nt^0 
=  1.68  -  0*i.  The  measured  and  calculated  data  from  500-800 
nm  at  all  angles  for  the  thiophene  film  are  plotted  in 
Appendix  I. 
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VI .  Recommendations 


Continuation  of  this  work  is  mandatory  to  show  the 
validity  of  the  models  chosen.  In  addition,  new  parameters 
may  be  added  or  new  approaches  tried  to  improve  fitting  the 
data. 

First,  the  number  of  data  sets  taken  at  different  angles 
of  incidence  should  be  increased.  Five  sets  of  tan  f  data 
only  fit  the  simplest  model  well.  Many  more  data  sets  are 
needed  for  more  complex  models.  Remember,  one  can  never  have 
too  much  data. 

In  order  to  improve  the  quality  of  data  taken,  long  range 
measurement  trends  must  be  observed.  This  will  give  an 
indication  on  how  far  to  trust  the  data,  as  well  as  giving  the 
user  insight  into  machine  performance.  Also,  to  use  more  of 
the  data  acquired,  the  compensator  assembly  that  came  with  the 
machine  could  be  used.  This  would  allow  the  use  of  the  cos  A 
data  >  I  0.7071  .  Another  technique  to  use  more  of  the  measured 
data  would  be  to  mount  the  films  on  an  absorbing  substrate 
(say  silicon) . 

At  some  point,  the  models  used  in  this  thesis  will  not  be 
adequate.  New  models  will  have  to  be  developed,  probably  to 
include  more  parameters.  Two  additional  parameter  sets 
immediately  come  to  mind: 

1)  Accounting  for  surface  roughness  of  each  film  layer. 

Much  work  has  been  done  in  assuming  an  Effective  Medium 
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Approximation  (EMA)  for  isotropic.  film(s)/substrates.  No 
one  has  developed  the  formulae  for  an  EMA  in  the 
anisotropic  case. 

2)  Accounting  for  spatial  dispersion  (inhomogeneity)  in 
the  f ilm(s)/substrate.  Again,  the  models  herein  assume 
homogeneity.  Once  spatial  dispersion  is  to  be  accounted 
for,  the  derivation  must  either  be  modified  or  a  new 
derivation  used. 

Quality  of  data  will  affect  both  of  the  above  mentioned 
changes.  I  predict  that  until  the  quality  of  data  typically 
produced  is  improved  by  at  least  an  order  of  magnitude, 
neither  of  the  effects  listed  above  will  be  accurately 
measured  for  anisotropic  films. 

Finally,  a  real  computer  programmer  should  take  these 
models  and  make  them  "user  friendly”.  A  shell  driving  all  of 
the  models  could  easily  be  written  by  such  a  person.  A  look 
at  the  code  written  for  this  thesis  shows  this  author  is  not 
such  a  person. 
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Appendices  A-D 


The  following  four  appendices  contain  the  four 
spectroscopic  ellipsometry  data  sets  used  in  this  thesis.  For 
each  angle  of  incidence,  tan  T  and  cos  A  data  was  taken  every 
5  nm  from  300-800  nm,  for  a  total  of  202  data  points. 
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tar  psi  cos  delta 


Appendix  A:  ITO  Spectroscopic  Ellipsometrv  Data 


Figure  7  ITO  SE  Data  at  60* 
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tan  psi  cos  delto 
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tan  psi  cos  delta 


Figure  10  ITO  SE  Data  at  75* 


tan  psi  cos  delta 


Figure  11  ITO  SE  Data  at  80* 
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Appendix  B:  PBLG  5  Laver  Spectroscopic  Ellipsometrv  Data 


Figure  12  PBLG  5  Layer  Film  at  60’ 
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tan  psi  cos  delta 
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ton  pal  coa  delto 


lambda  (nm) 


Figure  14  PBLG  5  Layer  Film  at  70* 
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tan  psi  cos  delta 


Figure  15  PBLG  5  Layer  Film  at  75 


tan  psl  cos  delta 
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Cos  Delta  vs.  Wavelength  (5  layer  PBLG  on  ITO,  80  degrees) 


lambda  (nm) 


lambda  (nm) 


Figure  16  PBLG  5  Layer  Film  at  80’ 
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Appendix  C:  PBLG  8  Laver  Spectroscopic  El  1 ipsometry  Data 


lambda  (nm) 


lambda  (nm) 


Figure  17  PBLG  8  Layer  Film  at  60* 
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tan  psi  cos  delta 


Figure  19  PBLG  8  Layer  Film  at  70* 
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tan  psi  cos  delta 
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tan  psi  cos  delta 


lambda  (nm) 


Figure  21  PBLG  8  Layer  Film  at  80* 
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Tan  Pal  Co*  Delta 


Appendi 


Thiophene  Spectroscopic  El lipsometrv  Data 
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Cos  Dslta 


Ton  Ps!  Cos  Delta 


Figure  25  Thiophene  Film  at  75* 


Appendices  E.  G-I 


For  each  of  these  appendices,  both  experimental  and 
calculated  tan  T  and  cos  A  data  is  plotted  versus  wavelength. 
For  Appendices  E,  G  and  H  the  data  used  in  the  fit  was  tan  Y 
from  305-370  nm,  for  a  total  of  28  data  points  at  each  angle 
of  incidence.  For  Appendix  I,  the  data  used  was  tan  T  from 
500-800  nm,  for  a  total  of  61  data  points  at  each  angle  of 
incidence. 

At  the  bottom  of  each  of  these  appendices  are  the  values 
used  in  the  final  fit,  as  well  as  the  variance  of  that  fit. 
Even  though  only  tan  Y  data  were  used  in  the  fits,  and  the 
variances  were  calculated  with  only  the  tan  Y  data,  the 
calculated  and  experimental  cos  A  data  is  also  plotted  versus 
wavelength  for  each  data  set. 
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endix  E:  ITO  Fitted  Data 


Figure  27  ITO  Fitted  and  Experimental  Data  at  60*  (Dashed  Line  is 
Fit) 


Tf  =  279  ±  3  A 

(N{  dispersive  from  Njqj  =  2.38  -  0.12*i  to  Njjq  =  2.15  -  0.07*i) 


VARIANCE  =  0.0093 


Wavelength  (nm) 


Figure  28  ITO  Fitted  and  Experimental  Data  at  65*  (Fitted  Data  is 
Dashed) 

T{  =  279  ±  3  A 

(N{  dispersive  from  Njjj  =  2.38  -  0.12*i  to  Njjj  =  2.15  -  0.07*i) 
VARIANCE  =  0.0093 
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Wavelength  (nm) 


Figure  29  ITO  Fitted  and  Experimental  Data  at  70*  (Fitted  Data  is 
Dashed) 


Tt  =  279  ±  3  A 

(Nf  dispersive  from  N3J5  =  2.38  -  0.12*i  to  Njjj  =  2.15  -  0.07*i) 
VARIANCE  =  0.0093 
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Wavelength  (ran) 


Figure  30  ITO  Fitted  and  Experimental  Data  at  75*  (Fitted  Data  is 
Dashed) 


T{  =  279  ±  3  A 

(N{  dispersive  from  Njqj  =  2.38  -  0.12*i  to  Njjq  =  2.15  -  0.07*i) 


VARIANCE  =  0.0093 


Wavelength  (nm) 


Figure  31  ITO  Fitted  and  Experimental  Data  at  80*  (Fitted  Data  is 
Dashed) 


T{  =  279  ±  3  A 

(Nf  dispersive  from  Nj^  =  2.38  -  0.12*i  to  Njjj  =  2.15  -  0.07*i) 


VARIANCE  =  0.0093 


_ Ton  Psi  vs.  Wavelength 
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lambda  (nm) 

Figure  33  Fitted  and  Experimental  Data  for  5  Layer  PBLG  on  ITO  Coated 
Glass  at  60*  (Dashed  line  is  fitted  data) 


T{1  =  164  ±  2  A 

(N,!  =  1.65  -  0 . 20*i  and 

Njj  =  =  1.62  -  0.10*i  for  all  wavelengths) 

VARIANCE  =  0.0075 
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(ombda  (nm) 


lambda  (nm) 


Figure  34  Fitted  and  Experimental  Data  for  5  Layer  PBLG  on  ITO  Coated 
Glass  at  65*  (Dashed  line  is  fitted  data) 


T{1  =  164  ±  2  A 

(N,!  =  1.65  -  0 . 20*i  and 

Nll  s  Nil  =  1,62  “  0.10*i  for  all  wavelengths) 
VARIANCE  =  0.0075 
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Figure  35  Fitted  and  Experimental  Data  for  5  Layer  PBLG  on  ITO  Coated 
Glass  at  70*  (Dashed  line  is  fitted  data) 


T{1  =  164  ±  2  A 

(N,!  =  1.65  -  0 . 20*i  and 

Nji  =  Njj  =  1.62  -  0.10*i  for  all  wavelengths) 
VARIANCE  =  0.0075 
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Figure  37  Fitted  and  Experimental  Data  for  5  Layer  PBLG  on  ITO  Coated 
Glass  at  80'  (Dashed  line  is  fitted  data) 


Tu  =  164  ± 
(Nlt  r  1.65 

=  N«1  = 
VARIANCE  = 


2  A 

-  0.20*i  and 
1.62  -  0 . 10*i 
0.0075 


for  all  wavelengths) 
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Appendix 


PBLG  8  Layer  Pitted  Data 


Figure  38  Fitted  and  Experimental  Data  for  8  Layer  PBLG  on  ITO  Coated 
Glass  at  60*  (Dashed  line  is  fitted  data) 


Tu  =  216  ±  2  A 

(Nfl  *  1.65  -  0 . 18*i  and 

Nyl  =  N,!  =  1.62  -  0.09*i  for  all  wavelengths) 
VARIANCE  =  0.0056 
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Glass  at  70*  (Dashed  line  is  fitted  data) 


Tfl  =  216  ±  2  A 

(Nji  =  1.65  -  0.18*i  and 

Njl  =  Nil  =  -  0.09*i  for  all  wavelengths) 

VARIANCE  =  0.0056 
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Figure  41  Fitted  and  Experinental  Data  for  8  Layer  PBLG  on  I  TO  Coated 
Glass  at  75*  (Dashed  line  is  fitted  data) 

Tfl  =  216  ±  2  A 

(Nxl  =  1.65  -  0 . 18*i  and 

Nll  s  Nil  =  -  0.09*i  for  all  wavelengths) 

VARIANCE  =  0.0056 
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Figure  42  Fitted  and  Experimental  Data  for  8  Layer  PBLG  on  ITO  Coated 
Glass  at  80*  (Dashed  line  is  fitted  data) 


T{1  =  216  ±  2  A 

(Nxl  =  1.65  -  0 . 18*i  and 

Nyl  =  Nil  =  I*®*  “  0.09*i  for  all  wavelengths) 
VARIANCE  =  0.0056 
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Tf  =  3060  ±  7  A 

(Nf  =  1.68  -  0*i  for  all  wavelengths) 
VARIANCE  =  0.0075 
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Figure  44  Thiophene  Fitted  and  Experimental  Data  at  65*  (Fitted  Data 
is  Dashed) 


T{  =  3060  ±  7  k 

(Nf  =  1.68  -  0*i  for  all  wavelengths) 
VARIANCE  =  0.0075 
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Wavelength  (nm) 


Figure  45  Thiophene  Pitted  and  Experinental  Data  at  70*  (Fitted  Data 
is  Dashed) 


T{  =  3060  ±  7  A 

(Nf  =  1.68  -  0*i  for  all  wavelengths) 
VARIANCE  =  0.0075 
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Cos  Delta  vs.  Lambda 


Figure  46  Thiophene  Fitted  and  Experimental  Data  at  75’  (Fitted  Data 
is  Dashed) 


T{  =  3060  ±  7  A 

(Nf  =  1.68  -  0*i  for  all  wavelengths) 
VARIANCE  =  0.0075 
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Figure  47  Thiophene  Fitted  and  Experimental  Data  at  80*  (Fitted  Data 
is  Dashed) 


Tf  =  3060  ±  7  A 

(Nf  =  1.68  -  0*i  for  all  wavelengths) 
VARIANCE  =  0.0075 
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